home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / IMTQL2.z / IMTQL2
Text File  |  1996-03-14  |  3KB  |  67 lines

  1.  
  2.  
  3.  
  4. ____IIIIMMMMTTTTQQQQLLLL2222((((3333FFFF))))                                                        ____IIIIMMMMTTTTQQQQLLLL2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      IMTQL2, SIMTQL2  -  EISPACK routine.  This subroutine finds the
  10.      eigenvalues and eigenvectors of a SYMMETRIC TRIDIAGONAL matrix by the
  11.      implicit QL method.  The eigenvectors of a FULL SYMMETRIC matrix can also
  12.      be found if  TRED2  has been used to reduce this full matrix to
  13.      tridiagonal form.
  14.  
  15.  
  16. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  17.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  iiiimmmmttttqqqqllll2222((((nnnnmmmm,,,, nnnn,,,, dddd,,,, eeee,,,, zzzz,,,, iiiieeeerrrrrrrr))))
  18.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
  19.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,,nnnn))))
  20.  
  21.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssiiiimmmmttttqqqqllll2222((((nnnnmmmm,,,, nnnn,,,, dddd,,,, eeee,,,, zzzz,,,, iiiieeeerrrrrrrr))))
  22.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
  23.           rrrreeeeaaaallll             dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, zzzz((((nnnnmmmm,,,,nnnn))))
  24.  
  25.  
  26.  
  27. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  28.      On INPUT
  29.  
  30.      NNNNMMMM must be set to the row dimension of two-dimensional array parameters
  31.      as declared in the calling program dimension statement.
  32.  
  33.      NNNN is the order of the matrix.
  34.  
  35.      DDDD contains the diagonal elements of the input matrix.
  36.  
  37.      EEEE contains the subdiagonal elements of the input matrix in its last N-1
  38.      positions.  E(1) is arbitrary.
  39.  
  40.      ZZZZ contains the transformation matrix produced in the reduction by  TRED2,
  41.      if performed.  If the eigenvectors of the tridiagonal matrix are desired,
  42.      Z must contain the identity matrix.  On OUTPUT
  43.  
  44.      DDDD contains the eigenvalues in ASCENDING order.  If an error exit is made,
  45.      the eigenvalues are correct but UNORDERED for indices 1,2,...,IERR-1.
  46.  
  47.      EEEE has been destroyed.
  48.  
  49.      ZZZZ contains orthonormal eigenvectors of the symmetric tridiagonal (or
  50.      full) matrix.  If an error exit is made, Z contains the eigenvectors
  51.      associated with the stored eigenvalues.
  52.  
  53.      IIIIEEEERRRRRRRR is set to ZERO       for normal return, J          if the J-th
  54.      eigenvalue has not been
  55.         determined after 30 iterations.  Calls PYTHAG(A,B) for sqrt(A**2 +
  56.      B**2).  Questions and comments should be directed to B. S. Garbow,
  57.      APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.